Infant Galaxy Clusters at Low Redshifts? 
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^ , Abstract 

! The population and population composition of galaxies in galaxy clusters at 
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present reflect the mass of the clusters and the mass growth of the galaxies in the 
past. We investigate them for six clusters. We show that galaxies in massive clusters 
stopped growing at redshifts of ~ 4. Moreover, we find that some small galaxy clusters 
>■ ■ (groups) have too many massive galaxies for their apparent masses. One possibility 

^ I is that these groups are much more massive and in a phase just before virialization. 

^ I If this is the case, they should be called 'infant galaxy clusters' that will be matured 

O ! clusters in the dynamical time-scale (~ 10^ yrs). 

O ■ Key words: galaxies: clusters — galaxies: halos — large-scale structure of universes 
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Clusters of galaxies are the most massive objects in the universe. It is generally be- 
lieved that dark matter constitutes a large fraction of the mass in the universe, and galaxies 
and clusters form in dark matter-dominated halos (dark halos). In the standard picture of 
hierarchical structure formation of the universe, small objects are the first to form and these 
then amalgamate into progressively larger system (Blumenthal et al. 1984). This means that 
a dark halo containing a galaxy cluster has formed via merging and accretion of smaller dark 
halos containing galaxies or galaxy groups (small galaxy clusters). From now on, we do not 
discriminate a dark halo containing a cluster from the cluster unless otherwise mentioned; the 
same rule is applied for a group and a galaxy. 

The mass distribution function (MDF) of galaxies absorbed into but not disrupted in a 
cluster is defined by the number of the galaxies in a given mass range and shows the population 
composition of the galaxies with different masses. The MDF tells us the average time when 
the galaxies stopped growing as follows. The mass of galaxies increases with time as matter 
accretes to them through gravity. As a result, the number of massive galaxies increases and the 
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MDF evolves. However, as the potential well of the host cluster (or its progenitors) deepens 
as the cluster grows up, the released gravitational energy of the cluster allows the galaxies to 
move fast in the cluster. The large velocity of the galaxies prevents them from gathering more 
matter and the galaxies cannot grow anymore. Thus, the MDF freezes after that. The MDFs 
observed at present should depend on the time of the freeze. On the other hand, the integrated 
MDF or the galaxy population in a cluster depends on the host cluster mass; the more massive 
a host cluster is, the more massive galaxies are contained in it. 

We have devised an analytical model to calculate the MDFs of galaxies in clusters (Fujita 
et al. 2002). In the field of cosmological study, the Press-Schechter (PS) model has been used 
to obtain an average MDF of objects (galaxies and clusters of galaxies) in the universe (Press, 
Schechter 1974). The PS model treats the evolution of initial density fluctuations of the universe 
and predicts the number of collapsed objects in a give mass range. However, this model cannot 
be used to obtain the MDFs of galaxies in clusters, because it does not take account of two 
important effects unique to clusters. One is that even before a cluster forms, the precluster 
region should have a higher density than the average of the universe. The other is the spatial 
correlation of initial density fluctuations of the universe. The correlation is important when 
the fluctuations are close to each other; a precluster region is an example of the crowded 
fluctuations. The former effect has already been involved in the so-called extended Press- 
Schechter model (EPS) (Bower 1991; Bond et al. 1991; Lacey, Cole 1993). The model used 
here is a further extension of the EPS model by taking account of the latter effect (the spatial 
correlation), so we refer to this model as the SPS model. The details of the SPS model are 
described in Fujita et al. (2002). 

2. Comparison between Theoretical Predictions and Observations 

In principal, we can calculate MDFs using the SPS model and compare them with 
the observed ones to estimate redshifts when the observed MDFs almost froze (that is, the 
formation redshifts of galaxies; Zf). We note that Zf is the average redshift when the growth 
of galaxies are held by external force (that is, the growth of the host cluster). At z = Zf the 
host cluster had not necessarily been unified, and the host cluster might be divided into several 
progenitors. We can say that the redshift Zf corresponds to the time when the average mass 
of those progenitors became large enough to hold the growth of most of the galaxies that are 
included in the host cluster at z = 0. Because it is not easy to measure galaxy masses directly 
by a method like gravitational lensing, it would be difficult to derive MDFs of galaxies in 
clusters observationally. Thus, we convert the MDFs derived using the SPS model to velocity 
distribution functions (VDFs) of galaxies by an empirical relation between the internal velocity 
and mass of galaxies obtained by numerical simulations^ (Fig. 13 in Bullock et al. 2001): 

^ The original paper assumed that the relation between the mass and internal velocity of galaxies is time- 
independent (equation [40] in Fujita et al. 2002). In this paper, we consider the time dependence as equation 
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log f 77^5^ ) = 2.09 + 0.0404^/ + 0.29 
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where Mgai is the mass of the galaxy and h is the Hubble constant represented by Hq = 
lOO/i km s~^ Mpc~^. Equation (1) can also be applied to clusters. We compare the VDFs 
with the observed ones. Here, 'internal velocity' means rotation velocity (late type galaxies) or 
dispersion velocity (early type galaxies). More massive galaxies have larger internal velocities. 
Moreover, there is a benefit to use VDFs instead of MDFs. The tidal interaction among galaxies 
in a cluster and that between a galaxy and the host cluster somewhat affect the MDFs and 
VDFs after their evolutions froze (Klypin et al. 1999; Okamoto, Habe 1999; Ghigna et al. 2000). 
This makes it difficult to estimate the time of the freeze from the MDFs or VDFs because 
the SPS model does not treat the interactions. However, recent ultra-high resolution N-body 
simulations showed that the VDFs are much less affected by the tidal interactions than the 
MDFs (Klypin et al. 1999; Okamoto, Habe 1999; Ghigna et al. 2000). We have confirmed that 
the VDFs derived using the SPS model are consistent with those derived by several simulations 
(Fujita et al. 2002). 

To compare the theoretical model with observations, we use the catalogues of early 
(E+SO) and late (Sp) type galaxies in seven nearby clusters (Giovanelli et al. 1997a; Scodeggio 
et al. 1998b). The clusters are A1367, Coma (A1656), A2634, A262, Pegasus and Cancer 
clusters and the NGC 507 group, although we do not discuss the Cancer cluster because the 
number of galaxies is too small (less than 10 within the cluster virial radius). We do not 
discriminate the NGC 383 group from the NGC 507 group because they are very close to each 
other. Cluster membership criteria are described in section 8 of Giovanelli et al. (1997a). In 
the SPS model, we fix the masses of the clusters at those obtained from X-ray observations if 
available (Reiprich, Bohringer 2002). The X-ray mass estimation does not suffer from the small 
number of galaxies especially for small clusters. Reiprich and Bohringer (2002) assumed that the 
average density of a cluster is 200 times the critical density of the universe. This is appropriate 
for the Einstein de-Sitter universe but not for the low density universe we adopted {Qq = 0.3, 
A = 0.7, (jg = 1.0, and h = 0.7). In this paper, we assume that the average density of a cluster is 
100 times the critical density of the universe, and the cluster masses in Reiprich and Bohringer 
(2002) are modified to be consistent with the assumption. On the modification, we assume 
that the dark matter distribution in clusters is approximated by a power law distribution with 
an index of —2.4 (Navarro et al. 1997; Horner et al. 1999). For the Pegasus cluster, since the 
mass is not obtained by X-ray observations, we estimate it by a relation between the velocity 
dispersion of galaxies in a cluster and the mass of the cluster; we assume that the relation is also 
given by equation (1) for Zf = 0. The velocity dispersion of the Pegasus cluster is 393 km s~^ 
(Scodeggio et al. 1998b). The resultant cluster masses (Md) and viral radii (r^ir) are shown in 
table 1. Note that these mass estimations assume that clusters and groups have collapsed and 
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virialized. On the other hand, we leave the average redshift when galaxies stopped growing in 
mass (zf) to be a free parameter in the SPS model. In order to deal with the observed rotation 
velocities of late type galaxies and the observed dispersion velocities of early type galaxies at 
the same time, we multiply the dispersion velocities by 1.5 to obtain the effective rotation 
velocities of the early type galaxies, because the rotation velocity of a late type galaxies (frot) 
and the dispersion velocity of an early type galaxy (a) are respectively given by 

<t^^,^'^-^, (2) 

'"gal (-^ '"gal 

where G is the gravitational constant, a is the constant (~ 2 — 3), and rgai is the radius of the 
galaxy. Thus, we define internal velocities as fint = frot for late type galaxies and Vi^t = l-5cr 
for early type galaxies. We assume that each dark halo contains only one galaxy. Actually, 
N-body numerical simulations show that a dark halo contained in a cluster does not contain 
smaller dark halos (subhalos in a subhalo) at least for dark halos of galaxy scales (Fig. 4 in 
Springel et al. 2001). This is also confirmed in calculations including an N-body simulation and 
semi-analytic model of galaxy formation (Okamoto 2002, private communications). Mostly, a 
dark halo of galaxy scales has only one galaxy. One may think that this is because of a poor 
numerical resolution. However, the time scale of dynamical friction for a subhalo in a galactic 
scale dark halo is 

t,„, = 3^3 X io» f " f ^] ' ( ,) ( ) " . (3) 

V 3 / V60kpc; V220kms-V V2 X IQIOMq/ ^ ' 
where InA is the Coulomb logarithm, r^h is the radius of the galactic scale dark halo, Wint is the 
internal velocity of the galactic scale dark halo, and Mgdh is the mass of a subhalo in the galactic 
scale dark halo (Binney, Tremaine 1987). This means that for a galactic scale dark halo, most 
subhalos in it should have merged with the galaxy at the center of the dark halo within the 
time scale of cluster evolution (~ 10^'^ yr) except for very small subhalos that would be difficult 
to be observed. We do not discuss dark halos of group scales 5 x 10^^ ^0) in clusters unless 
otherwise mentioned; those dark halos may contain multiple galaxies. The number of galaxies 
included in those halos is small (see §3). 

3. Results 

In Figure 1, we present the predicted and observed VDFs for A1367, Coma, and A2634 
clusters. These clusters are relatively massive (^ 4 x 10^^ Mq; table 1). The VDFs are rep- 
resented by cumulative numbers, that is, the numbers of galaxies having internal velocities 
larger than a certain value. The internal velocities are normalized by the velocity dispersions 
of galaxies in the host clusters (Vq; table 1), which are converted from the cluster masses using 
equation (1) for Zf = for consistency with galaxy internal velocities. The observational data 
are complete for Vi^t ^ 200kms~^ (Giovanelh et al. 1997b; Scodeggio et al. 1997). With the 
observational data, we count the two-dimensional numbers of galaxies on the celestial sphere 
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Fig. 1. Galaxy VDFs for Zf ~Q (dotted lines). Zf ~ 1.0 (solid lines), and Zf — 4.0 or 5.0 (dashed lines) 
for A1367, Coma, and A2634 clusters. The VDFs obtained by observations are shown by circles. 

within the angles corresponding to the cluster virial radii (A'^2d; table 1). Then, we multiply the 
number by 0.7 to exclude projected galaxies within the angles but outside the cluster virial radii 
in the three-dimensional space. The number 0.7 is appropriate when the galaxy distribution 
in clusters is approximated by a power law distribution with an index of —2.4 and has a cut 
off at several times the virial radii. The index of the power low distribution is consistent with 
observations and the results of numerical simulations if galaxy distribution follows dark matter 
distribution in clusters (Navarro et al. 1997; Horner et al. 1999). Most of the galaxies left in 
our samples (^ 70 %), especially for massive clusters, are early type galaxies that are mainly 
found in the central regions of clusters. Figure 1 shows that the VDFs calculated using the 
SPS model can reproduce the observations if the galaxies in those clusters formed at redshifts 
of about 4. The deviation of the observational data from the predictions at small velocities 
is due to the incompleteness of the observations of small galaxies. For the Coma cluster, the 
largest velocity point also deviates from the line. This can be attributed to the dynamical fric- 
tion that forces massive galaxies to fall into the cluster center and to merge with a central cD 
galaxy (Fujita et al. 2002). The derived formation redshift of the cluster galaxies is comparable 
to the redshift of major star formation of early type galaxies in massive clusters (Kodama et 
al. 1998). Since most of the galaxies in our samples are early type galaxies, this explicitly shows 
that early type galaxies in massive clusters are very old and have not grown for a long time in 
terms of their masses as well as the stars in them. Mass growth of the early type galaxies (e.g. 
mergers between galaxies) may be essential for the star formation. We note that if we adopt 
the Einstein-de Sitter universe as cosmological parameters (f2o = 1, A = 0, as = 0.7, and h = 0.5), 
smaller formation redshifts {zj ^2) are preferable to be consistent with the observations. 

Figure 2 shows the VDFs for A262, the Pegasus cluster, and the NGC 507 group. These 
clusters have relatively small masses 1 x 10^^ Mq; table 1) and we here call them groups to 
discriminate them from the above massive clusters. For a given Mnt/K); the mass of a subhalo in 
the groups is smaller than that in the massive clusters. In the hierarchical clustering scenario, 
less massive halos form earlier. This means that for a given redshift {z > Zf), hierarchical 
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clustering proceeds further on the mass scale of group subhalos than on the mass scale of cluster 
subhalos. Thus, for a large Vint/Vo and Zf, the predicted cumulative numbers A^(> V^nt/Vo) of 
the groups are larger than those of the massive clusters (Figures 1 and 2). 

Figure 2 shows that contrary to A1367, Coma, and A2634, the observed VDFs of the 
groups do not match the SPS model predictions regardless of Zf] these groups have too many 
massive galaxies for their masses. For late type galaxies, Giovanelli et al. (1997b) show that the 
observational data are complete for Vi^t ^ 200 km s^^ (A262 and NGC 507) and ^ 130 km s^^ 
(Pegasus). For early type galaxies, although the completeness of the data is not commented in 
the catalogue, we confirm that the data are complete at least to the completeness limit velocity 
of the late type galaxies from the magnitude distribution. The excess of observed galaxies in 
Figure 2 suggests that the actual masses of these groups are much larger than their apparent 
masses estimated on the assumption that they are relaxed and virialized. For these groups, it 
may indicate that matter including galaxies has already been gathered through gravity, but it 
has not virialized and relaxed to be matured massive clusters. That is, they may be 'infant 
galaxy clusters' just before virialization, although it is one of possible solutions. Figure 3 
actually supports this interpretation. This figure shows that the theoretical predictions match 
the observations perfectly if the group masses are larger than their apparent masses; the masses 
are assumed to be 7 times (A262), 10 times (Pegasus), and 16 times (NGC 507) larger than 
their apparent masses. The group radii and velocity dispersions are increased with the group 
masses (table 1). Strictly speaking, if these groups have not fully relaxed, it is not exact to 
estimate the galaxy numbers within the virial radii by multiplying the observed two-dimensional 
numbers by 0.7 as is the case where groups are relaxed. However, at least this estimation and 
the comparison with the SPS model predictions should be appropriate as long as the index of 
the galaxy distribution is not much different from —2.4. In table 1, A''2d of each group increases 
as its Mci and r^i,. increase. If we assume that the galaxy density distribution is represented by 
r~^'^, the number of galaxies within a radius r has the relation of oc r°'^. Since the projection 
of outskirt galaxies does not make significant difference between A^ and A^2D; the increase of A^2D 
is roughly consistent with this prediction. We adopt the lines of 2;/ = to be compared with the 
observational data because these groups are growing at present and thus Zf is expected to be 
close to zero. In fact, while the lines of 2;/ ~ 1 can also be adjusted to match the observations, 
the lines of a higher redshift {zj ~ 4) cannot because the slopes of these lines are different from 
the slopes of the observational data. This actually indicates that the galaxies formed recently 
in these group regions. In Figure 3, the deviation of the observational data from the predictions 
at small and large velocities can be explained by the incompleteness of observations of small 
galaxies and dynamical friction, respectively. 

We make a comment on subhalos in a subhalo here. Our model does not take account 
of galaxies in larger subhalos of galaxy group scales. However, we can show that the number 
of those galaxies is small as follows. We take A1367 as an example. The expected number of 
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Fig. 2. Same as Figure 1 but for A262, the Pegasus cluster, and the NGC 507 group. 




Fig. 3. Same as Figure 2 but the group masses are increased so that the model predictions match the 
observations. 

subhalos with the scales of galaxy groups (say Vi^t ~ 400kms~^) in this cluster is only about 
one (Figure 1). On the other hand, the internal velocity Vi^t of the subhalo is close to that of 
the Pegasus cluster (Vq; table 1). As is shown by the lines in Figure 2, the expected number 
of galaxies in a (sub)halo of this mass is ~ 1/10 of the expected number of galaxies in A1367 
for a given \^nt- This means that the galaxies in larger subhalos do not much affect our results 
shown above. 

4. Discussion 

If the groups we investigated are infant galaxy clusters, we expect that they will relax and 
virialize within their dynamical time-scales (~ 10^ yrs). Considering the selection effect that 
groups having many galaxies in small regions can easily be observed, many of groups observed 
so far may be such infant galaxy clusters. In particular, compact galaxy groups such as the 
Hickson groups (Hickson 1982), which are small galaxy groups that are crowded with galaxies, 
are good candidates of infant galaxy clusters. Some of these groups may be collapsing and 
virializing at present. In fact, observations support the idea. First, the properties of compact 
groups (e.g. the X-ray luminosities or the spatial distribution of the galaxies) vary widely 
among the groups (Pildis et al. 1995; Ribeiro et al. 1998). This may reflect that we are seeing 
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different dynamical stages of cluster or group evolution. Second, HCG 62, which is one of the 
Hickson groups, has hard X-ray emission (Fukazawa et al. 2001). This may be the result of gas 
accretion toward the group and the particle acceleration at the shock waves or in turbulence 
formed in the rapidly accreting gas; the accelerated particles could produce the nonthermal 
hard X-ray emission via inverse Compton scattering of photons (Fujita, Sarazin 2001; Fujita 
et al. 2003). In the future, direct observations of matter distribution around groups through 
gravitational lensing will be particularly important to reveal the process of cluster formation. 
For example, diffuse density enhancement may be observed around some groups. It will also be 
important to compare the results of the SPS model with those of gravitational lensing, optical, 
and X-ray observations for many clusters and groups systematically to reveal the whole picture 
of cluster formation. 

Finally, we note that Evrard et al. (1996) used cosmological gasdynamic simulations 
to investigate the accuracy of cluster mass estimates based on X-ray observations; they found 
that X-ray mass estimates are very accurate and the errors are 14%-29%. However, recent 
X-ray observations have shown that the structure of clusters is much more complicated than 
the results of numerical simulations done by Evrard et al. (1996) and have shown that many 
clusters have not fully relaxed. Ikebe et al. (1996) found that the Fornax cluster has two 
distinct spatial scales of the gravitational potential. Recently, Chandra observatory found that 
clusters often have substructures with small X-ray temperatures called 'cold fronts' (Markevitch 
et al. 2000; Vikhlinin et al. 2001; Mazzotta et al. 2001). These show that clusters often have 
dark halos with the scales of galaxy groups, which is also shown by recent ultra- high resolution 
simulations (Okamoto, Habe 1999; Ghigna et al. 2000; Springel et al. 2001). For the three groups 
we investigated, the observations of X-ray temperature (or velocity dispersion of galaxies) might 
have been made only for possible central substructures, because it is often difficult to observe 
the X-ray temperatures (or the velocity dispersion of galaxies) of the relatively diffuse outer 
regions of the possible host clusters, especially when the host clusters are small. This could 
lead to the underestimation of the host cluster masses. The hard X-ray emission from HCG 62 
may come from the diffuse hotter gas in the outer region of the cluster. If this kind of emission 
is observed in groups with an excess of galaxies, it would be an evidence that the groups are 
infant clusters. 

5. Conclusions 

We study the velocity distribution functions (VDFs) of galaxies in clusters. We show 
that galaxies in some clusters stopped growing at a very early time of the universe or redshift 
of ~ 4, which is comparable to the redshift of major star formation in early type galaxies in 
clusters. Moreover, we find that some galaxy groups have too many massive galaxies for their 
apparent masses. One possibility is that these groups are in a transient phase and are much 
more massive than their apparent masses. If this is the case, they should be called 'infant galaxy 
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clusters' that will be matured clusters in the dynamical time-scale (~ 10^ yrs). Considering 
the selection effect for observations, many of galaxy groups observed so far may be such infant 
clusters. 

I thank C. L. Sarazin, H. Yahagi, M. Nagashima, T. Kodama, S. Nishiura, and T. 
Okamoto for discussions. Y. F. was supported in part by a Grant-in-Aid from the Ministry of 
Education, Science, Sports, and Culture of Japan (14740175). 
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Table 1. Parameters for Clusters. 



Cluster 


(1014 Mq) 


(Mpc) 


^0 

(kins~i) 


iV2D 


A1367 


4.8 


2.0 


988 


41 


Coma 


16 


3.0 


1410 


103 


A2634 


4.5 


2.0 


971 


56 


A262 


1.2 


1.6 


661 


38 


A262t (x7) 


8.4 


2.4 


1163 


55 


Pegasus 


0.20 


0.69 


393 


13 


Pegasus^ (xlO) 


2.0 


1.5 


766 


19 


NGC507 


0.54 


0.96 


525 


27 


NGC5071' (xl6) 


8.6 


2.4 


1173 


39 



f The cases where the masses are increased. 
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